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Abstract: One- loop amplitudes are to a large extent determined by their unitarity cuts 
in four dimensions. We show that the remaining rational terms can be obtained from 
the ultraviolet behaviour of the amplitude, and determine universal form factors for these 
rational parts by applying reduction techniques to the Feynman diagrammatic represen- 
tation of the amplitude. The method is valid for massless and massive internal particles. 
We illustrate this method by evaluating the rational terms of the one-loop amplitudes for 
99 H,-fj ^ 77, gg gg,jj ggg and 77 7777. 
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1. Introduction 



The upcoming LHC experiment provides a strong phenomenological motivation Q to 
develop formalisms which allow the successful evaluation of partonic multi-leg processes 
at one loop. Recently, many different methods have been proposed to deal with this 
highly complex task. Apart from algebraic reduction algorithms, which for multi-leg pro- 
cesses lead to a proliferation of terms and motivate numerical |^, ^ ^, |5|, 0] or semi- 
numerical 1^ |9|, H^, |ll], |l^, 14, 16 1 treatments at some stage, twistor-space inspired 
methods [|l^, |l^, 21, 22, 23, have stimulated a great deal of activity to complete 

the task - initiated already more than ten years ago 1 25, 26, ^] - to evaluate one-loop 



amplitudes elegantly from their unitarity cuts |28, 30 1. These algorithms are fully suc- 
cessful for supersymmetric amplitudes or special classes of helicity amplitudes, where the 
UV behaviour is tamed, but for general Standard Model amplitudes, it is difficult to obtain 
information on the so-called rational polynomials, which are induced by the ultraviolet be- 
haviour of Feynman integrals. This was analysed in great detail by Bern, Dixon, Dunbar 
and Kosower in [^6| , where criteria for 4-dimensional cut constructibility of one-loop am- 
plitudes were derived. We will refer to these criteria, called "uniqueness result" in [^], as 
the "BDDK-theorem" in the following. 
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The application of unitarity cuts to calculate the cut-constructible part of non-super- 
symmetric amplitudes has seen a lot of progress recently, and lead already to some remark- 



able results, as for example the cut constructible part of the six gluon amplitude |24|. This 
amplitude, including the rational parts, also has been evaluated numerically at a certain 
phase space point |31]. 



Very recently, progress also has been made to determine the remaining rational ambi- 



guities within the unitarity-based method, using the so-called "bootstrap approach" |37 



38, 39, Ho]. These ideas have lead to the successful determination of previously unknown 



multi-leg QCD amplitudes gl|, ||, |4|] . 

An alternative approach to the evaluation of the polynomial terms of helicity ampli- 
tudes in QCD is worked out in detail in p2| , |33| , |3^ , where the missing rational parts of 
the so far unknown six gluon helicity amplitudes have been given. The authors have used 



reduction formulae in Feynman parameter space, based on the work of |35, |3^. Their 
formalism is designed for gluon amplitudes and massless internal propagators. 

In this article, we show that the rational polynomials can be obtained in a general 
way, for massive as well as massless amplitudes and arbitrary external particles, in isola- 
tion from the cut-constructible parts, using Feynman diagrammatic reduction techniques. 
Our method is based on the tensor reduction formalism presented in [^], where explicit 
representations for tensor form factors have been derived. Projecting on the ultraviolet 
sensitive part of these form factors, we obtain the rational parts. The projection leads to 
a considerable simplification of the reduction cascade, which yields relatively compact ex- 
pressions. This defines a method which allows for the automated evaluation of the rational 
polynomials of arbitrary one-loop amplitudes^. 

In this sense, the two approaches - construction of an amplitude by unitarity cuts and 
using Feynman diagrams together with tensor reduction - can be considered as comple- 
mentary: for the polynomial part, the Feynman diagrammatic approach seems to be more 
straightforward, while for the remaining parts of the amplitude the unitarity-based method 
often leads to faster and more compact results. Such a combination of techniques already 



has been employed in |45] to obtain the one-loop amplitude for a Higgs boson plus four 
negative helicity gluons. Of course, for such a combined formalism to be universally appli- 
cable, the generalisation of the unitarity based methods to massive internal propagators is 
required. 

The paper is organised as follows: First we give a definition of rational or non-cut- 
constructible terms of general one-loop amplitudes in section 2. As an illustration, we 
calculate the polynomial terms of some 3, 4, 5- and 6-point amplitudes in section 3. Section 
4 contains our conclusions. In the appendix we provide formulae which are useful for the 
extraction of rational polynomials of IR divergent amplitudes. 

2. Rational terms of one-loop amplitudes 

In this section we propose a definition of the rational part of a one-loop amplitude. To 
^Work on a similar subject has also appeared very recently in [ ^ . 
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disentangle rational polynomials of infrared and ultraviolet origin, it is convenient to define 
the rational part with respect to the corresponding IR regulated amplitude. As an IR reg- 
ulator, we use off-shell momenta for the external legs of our tensor integrals or masses for 
internal fermion lines. In this way it is guaranteed that all poles in e will be of ultraviolet 
nature and only terms which are related to ultraviolet divergences lead to rational poly- 
nomials in an amplitude. It turns out that in the resulting expressions, the limits to the 
original kinematics are well defined. Amplitudes which contain dimensionally regulated 
infrared divergences could be used directly as a starting point as well, however in this case 
the two- and three-point functions have to be treated in a different way, as will be discussed 
in the appendix. 

2.1 Definition of rational parts of amplitudes 

Each A^-point amplitude T can be written as a linear combination of tensor Feynman 
integrals in n = 4 — 2e dimensions. Schematically this can be denoted as 

r = ^C^,...^^(n,K-,mfc}) /;'^i-^«({s,„m,}) . (2.1) 

For our purposes, F can be regarded as a complex valued function of the dimension n and 
the kinematical invariants Sij , ruk ■ The tensor integrals are defined in momentum space 
as^ 

where qj = k + rj are the propagator momenta and rj = pi + • ■ ■ + Pj are sums of external 
momenta. In |12] we have presented all the relevant formulae to perform a complete tensor 



reduction of a general A-point tensor integral of rank R < N. For A^ > 6 the reduction 
is purely algebraic, in the sense that these rank R A-point functions decay into a linear 
combination of (A — l)-point functions of rank R-1. For A < 5 all tensor form factors are 
evaluated in terms of an adequate basis. We do not repeat these formulae here, but give 
the essential definitions to keep the paper self-contained. Introducing Feynman parameters 



leads immediately to the following representation of the tensor integrals |46, 47 



[R/2] / .Nm 

ji- Jr-. 



r^^^-^-{{s.,,m,}) = {-If (-^) E [(9 

X/]^+''"(jl . . ■jR-2m ; {Sij,mk}) . (2.3) 

The objects ■ ■ ■ ,jR-2m]{sij,'rnk}) are scalar integrals in D = n + 2m {m = 

0, 1, 2, . . .) dimensions with Feynman parameters in the numerator: 

lNUi,---jRi{sij,'mk}) = 

N N / N \ 

(_l)JVr(A - f ) / n dz. 6{1 - ^ zi) zj, ...z,J-Y, SMZkZi/2 

1=1 1=1 \ k,l=l ^ 

Ski = in - nf - mf - ml . (2.4) 



^To make contact to eq.(2.1) of we note that we have set = . . . = = here for ease of 
notation. 
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This leads to 



r = Y,C{n,{jl},{s^,,mk})I%-'^"'ajl}■As^J,mk}) , (2.5) 
where the sum runs over all different integrals making up the amplitude. 

Note that the coefficients C(n, . . . ) depend on the way the numerators of the Feynman 
diagrams are evaluated, i.e. on the renormalisation scheme which defines the dimension 
of the Clifford algebra and the dimensionality of the internal and external particles (for 



a more detailed discussion see |4^, |5^, 51]). If chiral fermions are present or if one wants 
to define helicity amplitudes for partonic processes, the presence of 75 makes it necessary 
to distinguish 4-dimensional from (n — 4)-dimensional contributions of the Dirac algebra 
and n-dimensional vectors. We note that due to the dimension splitting |^2| of the loop 
momentum k = k + k, where k is 4-dimensional and k is {n — 4)-dimensional, a few 
integrals with fe^-terms in the numerator have to be known, which can be mapped to 
higher dimensional integrals |27, ^ 



(k ■ kr _ , T{a-e)n-4 ^„+2a 



iW2(fc2-M2)^ ' r(l-e) 2 
f d-k (k ■ krk^^k- _ +ir(a-e) n-4n + 2a 2+2„ 

y iW2 (fc2 - M2)^ - ^ r(l-e)^ 4 n ' ^^'^^ 



Therefore, these integrals contribute to the sum in eq. (| 

The rational part of an amplitude in general stems from two different sources: Firstly 
from a linear combination of d = 4 terms with the finite rational terms of Feynman param- 
eter integrals, secondly from (n — 4)-dimensional remnants of the Dirac algebra and the 
treatment of internal particles, which combine with UV pole parts of Feynman parameter 
integrals. Thus we define the rational part TZ of an amplitude F as'^ 

7^[^] = C(4, {ji}. {si,,mk})n[Il+''"^{{3i}-. {si,,mk])] (2.7) 
+ (n-4) C\4,{ji},{s,,,m^})V[P+^^{{jiy,{s,,,mu})] 

with 

C'(4, {i«},{sij,mfc}) = ^C{n,{ji},{sij,mk}) 



71=4 



In eq. (p.7|), V is the projector onto the pole part of the argument, i.e. the 1/e-term in the 
e-expansion. The action of TZ on Feynman parameter integrals is outlined in more detail 
below. It is the coefficient C which governs the renormalisation scheme dependence. 

In [12 1 we have given explicitly all the necessary formulae to reduce the scalar integrals 
with nontrivial numerators to n-dimensional one-, two- and three-point functions, (n+2)- 
dimensional three-and four-point functions and (n-|-4)-dimensional four-point functions. 



^We assume that F is either a genuinely UV finite or an UV renormalised amplitude. Note that coun- 
terterms can be expressed in terms of one- and two-point functions. 
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Schematically 

/^''"(U}) = E ^1 + E i2inji\ji,h) (2.8) 

+ E /?3a({iJ)/3"(l|jlbl,J2|jl,j2,i3) +E /53b({ia)/3"+'(l|il) 

+E /?4a({ij)/r'(iijiiii,i2iii,i2,i3) +E /546(j)/r'(iii)]} > 

where the arguments j2lii, J2, J^) denote integrals with up to three Feynman pa- 

rameters in the numerator, and are the genuine iV-point scalar integrals, which will 

be denoted simply by in the following. These integrals we call the "GOLEM integral 
representation'"^. It is defined by the property that inverse Gram determinants can be 
completely avoided by using such a representation. The coefficients Pk{{ji}) are polyno- 
mial in the kinematical variables {sjj,m|}, they do not depend on the dimensionality n. 
The GOLEM integral representation is a preferable starting point for a numerical evaluation 
of one-loop amplitudes. In algebraic approaches it is useful to reduce the GOLEM integrals 
further to a smaller integral basis which allows for an easy isolation of IR/UV divergencies. 
A convenient choice is to express each GOLEM integral by a linear combination of the scalar 
integrals If, I^, I^, 1^+^. AH necessary formulae can be found in ||l2[. It is only in this 
further reduction step that an n-dependence enters into the coefficients: 

^r'^ciij) = E ^i(^) ^" + E I2 + E ^3(^) + E ^r' • (2-9) 

The summation is understood over the different kinematically allowed 1,2,3- and 4-point 
functions, which are defined by all possible propagator pinches of the corresponding Ap- 
point function on the left-hand side. 



To completely define the operator TZ introduced above, we need to determine its action 
on a linear combination of 1,2,3 and 4-point functions of the form given above. Using again 
the rule 

n[c{n)lN] = c(4) n[lN] + (n - 4) c'(4)P[/,v] , (2.10) 

one sees that the rational part of an arbitrary integral of type is defined by 

the rational and the pole part of scalar integrals with trivial numerators. As UV poles can 
only occur for D/2 + m > N , all D=4 three-point and D=6 four-point functions are UV 
finite, i.e. their pole parts are zero: 

V[l!] = , nil] = . (2.11) 

Further, the BDDK ||2^ theorem tells us that all polynomial terms of such integrals, e.g. 
terms ~ tt^/G, are fully reconstructible by considering 4-dimensional cuts. This suggests 
to define the rational parts of these integrals to be zero: 

7^[/|] = , 7^[/|] = . (2.12) 
* GOLEM stands for "General One-Loop Evaluator of Matrix elements" ||l6[. 
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A detailed description of how these functions are uniquely reconstructible from the asymp- 
totic logarithmic behaviour of an expression can be found in |25, 26 1. 

We define the pole- and rational parts of two-point functions for s 7^ by considering 
the following massive representation 

1 



I2 {3,1711,1712) = / dxlog{—sx{l — x) -I- xml + (1 



x)ml) 



as 



V[I^] = i , nr^] = . (2.13) 

Note that in the case rrii = 0, another natural definition could be 7^[/2 ] = 2, but these 
rational terms are directly related to the cut-constructible logarithmic terms. 

In the case s = the two-point function degenerates to combinations of one-point 
functions. The different cases are related to one-point functions in the following way: 



mf — rrin V / 



/2"(0,0,m2) ^ 2 ^ ^2 ^ 2) _ (2.14) 

n — 2 



As up to 0(e) 



I"(m^) = ^1 - log(m^) , 

we define the the pole and rational part of the one-point function to match the non- 
logarithmic term, i.e. 

2 

p[/«(m2)] = , n[I^{m'^)] = . (2.15) 

Note that the rational parts of 12(0, 0, m^) and /2 (0; mPjrn^) turn out to be different using 
this definition, on the other hand they still respect relation (|2.14 ). 

After having defined the pole and rational parts of the scalar integrals I^'^'^ , -^3 , I2 



, the rational part of an amplitude, Eq. (2.7), is fully determined. For a renormalised, 
IR-finite amplitude, we are now in the position to define the cut-constructible part of the 
amplitude T indirectly by 

c[r] = {1 - n)[r] . (2.16) 

In the general case one may use the definition 

c[r] = lim(l-7^)[^IR 

regulated ] • (2.17) 

IR 

By an "IR regulated" version of the amplitude we mean here a representation where in- 
ternal propagator masses or virtualities of external particles which are zero in the original 
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amplitude are treated as non-zero, i.e. nij 7^ or sj = pj ■ pj 7^ 0, in tensor form factors. 
This renders the amphtude IR finite in the hmit n ^ 4 and the extraction of pole and 
rational parts is as in the IR finite case, if we demand limjR Fir regulated = T, where limiR 
denotes the limits nij — > and/or sj 0. In doing so it is of course crucial to avoid 
the expansion in (n — 4), which does not commute with the limits rrij 0,Sj — > 0. The 
IR-limits ^ or s^- ^ do not commute in general with the extraction of pole and 



rational terms. An example will be discussed below in subsection 3.3. In the appendix we 
provide explicit expressions for pole and rational parts of IR divergent tensor form factors, 
which can be used to define rational terms differently, without using this IR regularisation 
procedure. 

We have checked that in the case of massless internal particles our definition of cut- 
constructibility leads to identical form factor expressions as the ones presented in 



which were used to confirm known results for five- and six-gluon amplitudes derived with 
the methods of ||2^ . We stress that it is not necessary to have an off-shell representation 
of the full amplitude; the IR-regulated representations of the needed tensor form factors 
are sufficient to define the rational and pole parts we are looking for. 
For later use, we also define the following operator: 

U[I] = {V + n)[I] . (2.18) 



Note that the unitarity based methods ^ are not yet generalised to the case 



of massive loop integrals. Our definition can easily be adapted to modified definitions of 
cut-constructibility for massive one-loop amplitudes by redefining rules (p.ll|) -( ^T5|) , once 
such a formalism is developed. 

For practical purposes, it is more convenient to produce purely rational expressions for 
tensor integrals directly once and for all, instead of reducing to scalar integrals first and 
then apply the operators V and TZ, because this avoids an explosion of terms at intermediate 
stages of the calculation. Therefore, we give a list of the pole- and rational parts of the 
form factors for tensor integrals for = 2, 3 and 4, which were needed in the applications 
below. As will be explained below, their knowledge is sufficient for the determination of 
the rational parts of any A^-point amplitude, including amplitudes with massive internal 
particles. 

2.2 Rational parts of 2-point form factors 

Higher dimensional 2-point functions can be reduced to (4 — 2e)-dimensional ones by ap- 
plying scalar integral reduction formulae |35, 36, 47, |32| . 



X {-s + m\- ml)Ii{ml) + {-s -m\ + ml)Ii{m\) + A(s, ml, ml) 12(3, m\,ml) 
A(x, y, z) =x^ + y'^ + - 2{xy + yz + xz) . (2.19) 
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Applying lA to this formula and using rules ( ^.13 ),( 2.15| ) yields, for the special case m\ = 

(2.20) 



2 2 
777-2 ~ 



m 



1 + 1 

e 



6 Ve 3/ 

The rational parts of the 2-point form factors can be read off directly from the formulae in 
appendix A of [12|. The form factors are defined by 

For equal internal masses one obtains 

j2,2 



l2,2/ 





1 






2^ 






s 






12 




1 










+ 


3 





2\ 

+ 3 +T 



+ 1 



(2.21) 



Formulas for different internal masses are obtained analogously, but are not listed here as 
they are rather lengthy and because they are not needed in the following. 

2.3 Rational parts of 3-point form factors 

Again, the rational parts can be obtained by applying U = TZ + V as defined above to the 
reduction formulae in section 5.1 of [^]. In order to give compact formulae for the rational 
part of the three point form factors, we will introduce the following matrix: 



Hi 



h h 



B 



SfX ,hJ2e {1,2,3} 



(2.22) 



where 6/ = Sjj} and B = Ylk=i ^k- Note that internal propagator masses are present 

through the matrix S defined in eq. ( p.4[ ). Applying momentum conservation = pi + 
P2 + P3 = and defining Gik = 2 ri ■ iov l,k € {1,2}, it is easy to see that the third 
minor of H is just G^^: 



2 S3 



S2 - Si - S3 



(2.23) 



A(S1,S2,S3) I S2 - Si - S3 2 Si / 

with Sj = pj -pj and A(si, S2, S3) as defined in eq. (|2.19|) . Then one can define the quantity 



Hhh 



1 + 6, 



+ 



h h 



+ 5Hi,i,^ +1^3 + 2^3 

JD 



.(2.24) 



Here the ratios bj/B, where 



bi _ -S3 (si + S2 - S3) + m| (si - S2 - S3) + 2 S3 mj - m| (si - S2 + S3) 



B 



A(S1,S2,S3) 
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and B,h^/ B are obtained by cyclic permutations of indices, are well defined in the 
limit of massless internal propagators and light-like on-shell kinematics, as they behave 
like 1/A(si, S2, S3), which is well defined as long as at least one Sj is non-zero. Thus H 
and V are both well defined for all relevant kinematical cases, i.e. m| — > and/or Sj — > 0. 
Note however that these limits are not always equal to the finite parts of the corresponding 
integrals in the case where one or two of the Sj are vanishing, as on-shell limits and e- 
expansion do not commute in general. This issue will be treated in detail in appendix 

0- 

The tensor form factors are defined by 

2 

I3^^{ri,r2,rs = 0, mi, 7712,7723) = ^ Af^^ r'^^ 

ii=i 



2 



I«'A'1M2 ^ (^^,^2,r3 = 0,7771,7772 , 7773) = i?^'^ 9^'^' + ^ Af^^ r^W^^ 

iij2=i 



-^3 ^ ^ ^ = (n,r2,r3 = 0, mi, 1712, m^) = 

2 2 

^53,3^^^,^,^7^3 + 2perms.) + ^ ^S^js^n^'n'-js ^ (2-25) 

jl = l ilj2,J3 = l 

where we used momentum conservation to have = 0. Using the equations in section 5.1 
of ref. [12 1, one gets the following pole- and rational parts of the form factors: 



mi''] 


= 


U[B''^] 




mf;l] 




U[Bf'^] 


= --(-- 

12 Ve 




= -^khh ■ 



(2.26) 



2 k 
- + — 

3 ^ 



(2.27) 



The rational part of the rank one tensor integral is identically zero. In the massless case 
these formulae are identical to the ones derived for massless internal particles in [^ ]. 

The rational parts of the (77 + 2)-dimensional three-point functions are implicitly de- 
fined, e.g. for the case of all masses equal one finds 

U[I^+Hsi,S2,S3,m\m\m^)] =U[-2B^^^] = -^(^-^+l^ . (2.28) 
2.4 Rational parts of 4-point form factors 



From the preceding subsection and the way the form factors have been computed in ref. 112], 



it is clear that only the rank 3 and rank 4 form factors can have a non-zero rational part. 
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The form factors are defined by |12| 



^ 54,3 (^,,,,^.3 + 2 perms.) + ^ A^5,,, r^^r^^r^^ (2.29) 

jl = l ilj2j3 = l 

jn, ^llfJ,2^l3^^^ 
-'4 

3 



(ri, r2, r3, = 0, mi, m2, ms, 771,4) = C^''^{g^^^'^ g^'-^^"^ + 2 perms.) 



+ E (^^^'^^-^r-^ + 5 perms.) + ^"nnn ^n^n'n^^ -(^-^O) 



For rank 3, we get 



4,3i 



(2.31) 



and for rank 4: 



'l'2'3'4-l 



1 



24:6^ 72 
1 



12B ^ 



f^'^^ihjh; h, U) + /"^'^(^i, ^3; ^2, ^4) + /"^'^(^i, ^4; ^3, ^2) 
+/^''^(^2, ^3; ^1, ^4) + f^'^ih, h; h, h) + f^'^^ih, h; h-,h) 



(2.32) 
(2.33) 



(2.34) 



/'^''^(^i, ^2; ^3, k 



5^''^(^i; ^2, ^3, U) 



Ojl 



12 B 



h bj + 2 hi Hj (2 + 2 ^2 j 



{j} 
I3 U 



Hi 3^1 
Oif j = / 



In these formulae S is the index set labehing the internal propagators of the four-point 
function. H^^^ and V^^^ are related to the 3-point kinematics obtained when omitting 



propagator j from this set |12]. For completeness we also list 



1 5 

6e 18 



(2.35) 
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2.5 Rational parts of 5-point form factors 

With tlie results given above and the explicit representations of the 5-point form factors 



given in |12], it is manifest that up to rank 3, all 5-point functions have vanishing rational 
terms: 

U[I^] = U[ll^'] = Ull^'^'f"'] = ^ . (2.36) 

The rational terms of the form factors for the rank 4 and rank 5 tensor integrals can 



directly be obtained from the explicit formulae given in section 6 of |12], which manifestly 
show that no integrals apart from the ones given in the previous subsections appear in 
the reduction. For > 6, the reduction is purely algebraic anyway, i.e. involves only 
kinematic matrices until five-point functions are reached. Therefore, the formulae given in 
the previous subsections are sufficient to calculate the rational parts of arbitrary A^-point 
amplitudes. 

As was shown in |12], A^-point functions up to rank A^-2 are algebraically reducible 
to rank 3 five-point functions, therefore all the corresponding rational terms are zero. 



We have thus re-derived, within our formalism, a well known result of |26|: all A^-point 
amplitudes which contain, in a convenient gauge, at most rank A^-2 tensor integrals are 
cut-constructible. 



3. Applications 

In the following, five examples will be presented to illustrate our approach. The first is 
Higgs production by gluon fusion, the second scattering of light-by-light. Then we consider 
the 4-gluon amplitude which is IR and UV divergent. Finally we discuss a pentagon and 
a hexagon amplitude, ggg^^ —> and 777777 0. 

3.1 Example 1: Higgs production by gluon fusion 

Higgs production by gluon fusion is mediated by massive quark loops. In the Standard 
Model, the top quark provides the leading contribution. Up to a trivial colour structure 
~ 5ab the amplitude is given by 

M = - — [ ^"^ tr(^i(^i-Fmf)(^2+?TT-t)^2(^ + ^t)) 1 n 

V (47r)"/2 J i W2 [ql - m1){ql - mDik"^ - m?) ^ ' ' 

Here Qj = k -\- rj with ri = —pi and r2 = P2, where pj are the light-like momenta of 
the gluons with polarisation vectors £j. Apart from m| the only non- vanishing Lorentz 
invariant variable is 2pi • P2 = s = m?^ . Working out the trace leads to 

M ~ /g (ri,r2,0,m?,m?,m?) {-ei • £2 + 2 ei • P2 £2 • Pi + 2mlei ■ £2) 

+l!^^''{ri,r2,0,mf,mf,mf) {8ei^e2u - 2ei • £25^*^^) (3.2) 

The scalar integral /g does not contribute to the rational part of the amplitude, while the 
rank 2 tensor integral does. The pole and rational parts of the form factors for I^^'^ defined 
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above turn out to be 



3,21 _ 1 + e 



"1'^ ■ 4. 

U[A'^i\ = -U[AS\ = U[aS\ = . (3.3) 

The decomposition of the Feynman diagram in terms of partly UV divergent tensor integrals 
made it necessary to work in dimensional regularisation, although the amplitude is finite. 
The rational parts of the amplitude are a result of products of UV 1/e poles and order e 
remnants from the n-dimensional gamma algebra. Note that the rational part of the tensor 
functions is only a small part of the whole tensor integral. Applying the rules of section 2, 
the rational part is found to be 

n[M] = ^ ^ . (3.4) 

TT V S 

Here !F'^'^ = pjEj — pjEj is the abelian part of the gluon field strength tensor. 
3.2 Example 2: Scattering of light-by-light 

In QED, scattering of light-by-light is mediated by a closed electron loop. It is well known 
that the six box topologies making up the amplitude are related, such that it is sufficient 
to evaluate only one diagram with a given ordering of the external photons. The others 
can be obtained by all non-cyclic permutations of the photon momentum and polarisation 
vectors. 

-^ = 7i^ Z2 ^(^1' ^2,^3,^74) (3.5) 

In addition, due to charge invariance, one has ^(1,2,3,4) = t/(l,4, 3,2), ^(1,3,4,2) = 
2, 4, 3) and 4, 2, 3) = Q{1, 3, 2, 4). We evaluate now the rational polynomial of A4 
using the polynomial tensor coefficients. The evaluation of the diagram 

g(l 2 3 4) = - / ^"^ tr(^i(^i+ me)^2(^2+ me)^3(^3+ me)^4(^ + me)) 

involves four-point tensor integrals up to rank four. These are in general complicated 
functions, but the rational polynomials are actually very simple. Note that we are only 
interested in the massless case rue 0. We only keep the electron mass as an infrared 
cutoff for the moment. Otherwise the massless on-shell 2-point functions, which are zero in 
dimensional regularisation, would spoil a clear separation of IR and UV problems. We give 
now the complete list of the polynomial part of the four-point tensor coefficients. For rank 
zero, one and two, no rational terms are present. The rank 3 and rank 4 tensor coefficients 
are defined in eq. ( |2.29| ) above. The rational polynomials of the tensor coefficients depend 
on the external kinematics. In our case, where pj ■ pj = for j = 1, 2, 3, 4, we find in the 
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limit nip — > 



4,3i 



U[B, 



^[^333] 
^[^33] 

u — s 
2stu ^ 



u — t 

2stu 
1 



2su 



-U\A 



4,3 1 
I23J 



2tu 



(3.7) 



the remaining ones are defined by symmetry under exchange of the lower indices. The non 
cut-constructible parts of the rank four form factors are 



'4,4] _ J_l I A 

24 e 72 



,4,4i 



^[^1122 

^[^1133 
^[^1222 
^[^1223 



■ U\A 



4,4 1 
3333 J 



1 



= U[Blt] 

1 1 
— + 



UiBt^'] = -U[B',i] 



st su 2u^ 



4,4 1 _ 



■ ^[^2333 

- ^[^2233] 



1 



2su 
1 



1 

2^ 
1 



2st 2su 
1 1 

+ 



+ 



2u' 



6st 
1 

+ 



2m2 

1 



1 



6st Qsu 2v? 



1 



3st 



■ + 



1 



" ^[^2223] ~ 



3 s n ' 2 
1 



1 



1 



2st 2su 2v? 



1 



+ 



1 



+ 



1 



6si Qsu 2v? 



12 u 



(3.8) 



the remaining ones are defined by symmetry. Evaluation of the rational part of eq. (|3.6D is 
now straightforward. Using these form factors we find for the rational part of the sum of 
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all graphs the manifestly gauge invariant result 

^[ X] ^(^^1' '^2,0-3, 0-4)] 

(TG54/Z4 

8 tr(^i^2) tr(^3^4) ^ 8 tr(^l^tr(^b^ ^ 8 tr(^i^4) tr(^2^3) 



3,s s 3 u u 3 t t 

64 ti{TiJ^2) P4 • -^3 • Pi Pa • -^4 • Pi 64 tr(J^iJ'^3) pi ■ ■ ps P3 ■ ■ pi 



3 s 
64 tr(Jf^ijr4) 




stu 

■P3 PA 


•^3 


3 n 

Pi 64 tr(:r2jr3) p2 




stu 

P3 P3 


J^4 -Pl 


3 t 

64 tr(^2^4) P2 




stu 

■P3 P4 


^3 


3 t 
pi 64 tr(J^3j^4) p2 




stu 

P3 Pi 


•^2 •P3 



3 u stu 3,s stu 

_^-^Q24 P2 • -^1 • P3 Pi • -^2 • P3 P4 • -^3 ■ Pi P3 • -^4 • Pi g-^ 

where J^j'^ = Pj — Pj£j is the electromagnetic field strength tensor. On the amplitude 
level the UV pole cancels and thus TZ = IA. The result simplifies further if one specialises to 
helicity amplitudes using spinor helicity methods |53|. Due to parity invariance and Bose 

symmetry only the helicities + + ++, + + H — and + H have to be considered. 

The rational polynomials of the three helicity amplitudes are then given by 

„2.+ .+ _ o„,2 [12] [34] 



n[M++++]= 8a^e+ -£+4 •£+ = 8 a' 



(12) (34) 

7^[7\/^+++"] = Sa'' et ■ et et ■ piej -pi — = -8a 



2.+ 2. _ , [12] (14) [13] 

■^2^3 -Pi^. -PI- --8a ^^2) [14] (13) 

7^[X++-] = -8 a^ et ■ e+ £3 ■ = -8 a^ . (3.10) 

We note that up to a phase the rational parts of the different amplitudes are the same. For 
the first two helicity configurations, the rational terms are the full result. For completeness 



we also quote the full result of the + H case [54, 55 1: 

(12) [34] 



-B«Mi+^io.m+^^k^m+^^ii Mi. (3.11) 



3.3 Example 3: Gluon-gluon scattering 

To apply our formalism to an IR divergent amplitude we have chosen the 4-gluon amplitude 

as an example. In the following we consider the helicity amplitudes + + ++ and + H . 

Using spinor helicity methods one can replace the polarisation vectors by polynomial terms 
in the external momenta times a global phase [tr^(. . . ) = tr(. . . )/2 it tr(75 . . . )/2]: 

^+..^+..^+.3,+ .4 ^ ||[l|_Ltr-(^n'^^^2T^i)tr-(^37''^^^y'^3) 



-/.3.-/.4 _ [21] (34) 1 



(12) [43] 2s2 



tr- ( l^rf^' ^27^' ) tr+ ( I^Tf^^ /^^y^^ ) . (3.12) 



As we are using the 't Hooft-Veltman scheme for the computation, these expressions im- 
plicitly project onto the 4-dimensional part of n-dimensional objects in loop diagrams. At 



- 14 - 



tree level one finds 



with the colour objects defined through a trace of colour matrices in the adjoint represen- 
tation r„^^ = -i/^'^'' 

We have computed the next-to-leading order contributions in two ways to outline the 
extraction of the rational terms for the IR divergent and IR regulated case. As the + + ++ 
case is IR and UV finite, both evaluations give the same result for this helicity configuration: 

•^NLO ~ ''(12) (34)^^° 



-^NLO — 3'.^ad + -'ad + -'ad ) ■ (^CS.iOj 



The amplitude is defined by rational parts only. In the + H case we again extract an 

overall factor: 

- i^2 [21](34) 
•^NLO ~ (12) [43]^'"'-' 

i_ri- + Ax- + (^ + i)ri-). {3.16, 

In our definition we call the last line the rational part of the amplitude. Note that other 
polynomial contributions emerge after expanding the scalar integrals in e. We have checked 
that our result is — up to trivial factors stemming from different conventions — identical 



to the amplitude representations provided in ||57 



If we compute the same amplitude using off-shell values for the external momenta, i.e. 
kf = k2 = k'^ = = m?, some scalar functions change to their off-shell counterparts, 
e.g. 13(5) {s , m? , rn^) , I2~^'^{u,t) I2~^'^{u,t,m'^,m'^,m'^,m'^), which are all IR finite 
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now. The coefficients of the IR regulated basis integrals and the constant part are up to 
terms of order O(m^) identical to the coefficients of the original basis integrals. The only 
difference stems from the no-scale 2-point functions which are replaced by their off-shell 
counterpart, /2 (0) ~^ l2{rn?)- This leads to an additional term which is proportional to 
the leading order amplitude: 



+ -e + 0{m^) ] I^m^) . (3.17) 



32 4 
"3 

This additional term vanishes in the on-shell limit, as ^2 (0) = in dimensional regularisa- 
tion. When taking the on-shell limit it has to be put to zero before expanding in e, i.e. before 
the pole/rational part is extracted, otherwise the limit does not exist. This result defines 
the IR regulated version of the 4-gluon amplitude. The IR-limit limiR {Tin regulated) = T 
is smooth as long as scalar integrals are not expanded in e. For 2-point functions with 
limiR = the rational/pole parts must not be isolated from the result before the limit 
is taken. In practical calculations it is very easy to take care of these terms separately. This 
reasoning shows that the rational part of an amplitude as defined above is not affected by 
IR divergences. 

3.4 The ^jggg —>■ amplitude 



Following 1 56] the helicity amplitudes can be written as 



n2„3 

^AiA2A3A4A5 _ _fJ^yC3C4C5^AiA2A3A4A5 -j^g\ 



Due to Furry's theorem only one colour structure ~ f^^c^cs exists. The diagrammatic 
structure implies that this amplitude is not cut constructible, as rank four and five 5-point 
functions and rank three and four 4-point functions are present. Six independent helicity 
amplitudes exist. For three of them the cut-constructible part is identically zero, they are 
given by the rational part only: 

2 S34S45S35 

Here J^j"^ = pjEj — pjSj is the abelian part of the gluon field strength tensor. 
++++ _ tr(.^+^+)tr(.F+.F5^ 



A 

with the coefficient 
Further 



„2 2 
*23*45 



C-++++ p2 ■ • :P4 - (4 ^ 5) (3.20) 



= _£i5£i2 _ £15 _^ £23 _ £15 ^3 
■S24S35 535 S24 S34 



„2 2 
*12'^34 



pi . T7 . p3 - (3 ^ 4) (3.22 
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with the coefficient 

_ _£45£l3£l4 _ ■513^45 _^ ^Is _ ■^12 + '^45 ~ ■^12.S45 _^ glSSlS _^ glS - ■S34 
S35'Sl5S24 S15S35 S15S24 S35S15 S23S45 S23 

^34^45 _^ gl5 - ^25 _ ^23 + ^35 _ ^23^25 _^ ^34 + ^12 
S23'S15 S45 Sl3 S13'S45 ^15 

The results for the hehcity amphtudes A A^^^ and A contain also con- 

tributions from the cuts. The full result for these amplitudes can be found in [^]. We have 
verified that we get the same result for the rational terms, called A^ A^^^ and 

Ai^~^^ there, using our algebraic implementation of the rational polynomials defined by 
our formalism. 

3.5 The 6-photon amplitude 

Due to Bose symmetry and parity invariance, only four independent helicity amplitudes 
have to be evaluated, out of which two, the "all plus" and the "one minus" amplitudes 



identically vanish [59|. We evaluated the 6-photon amplitudes along the lines of the 4- 
photon case, with the difference that after taking the trace, all reducible scalar products 
in the numerator were cancelled directly. We have verified that the rational parts of 
and evaluate to zero using our formalism. For the non-vanishing 

amplitudes, we find 

7e[7\/[++++— ] = (3.24) 
■R[M++^—] = , (3.25) 



where eq. ( p. 24 ) is already known from the analytic result [^0|. The results have been 



obtained by two independent calculations, one based on the approach outlined in section 
2, the other one based on IR divergent form factors (see appendix]^). We observe that for 
the evaluation of the rational parts, form factors for at most rank 4 four-point functions 
were needed. Kinematically they are of the same complexity as the ones used for the recent 
evaluation of the rational parts of the six gluon amplitude [^] . 

4. Conclusions 

We have presented a formalism to evaluate the rational polynomials of arbitrary one-loop 
A^-point amplitudes. It is based on a tensor form factor representation derived in The 
definition of rational parts of these form factors induces a definition of the rational part of 
the amplitude. 

To disentangle contributions form UV and IR poles, we define ffist the rational poly- 
nomials with respect to IR-regulated amplitudes. We obtain compact expressions for the 
rational and pole terms of the tensor form factors which allow for an on-shell limit after- 
wards. In this way it is clear that the polynomial part of an amplitude origins only from 
the UV behaviour of the amplitude. This procedure defines rational polynomials of general 



one-loop amplitudes corresponding to the definition advocated in the literature ||2^, 32]. In 



addition, we give all the formulae to work with form factors which contain IR divergences 



-17- 



in the appendix. For IR finite amplitudes both approaches obviously lead to the same 
result. 

Both approaches were implemented in algebraic manipulation programs to allow for 
the fully automated evaluation of the rational parts of one-loop amplitudes starting from 
Feynman diagrams. The formalism has been applied to the evaluation of the amplitudes 
for gg — > iJ, 77 — >■ 77 , g'g' — > gg, 77 ggg and 77 ^nil- For the first four examples we 
recover the well-known results. For the six-photon amplitude we have proven by explicit 
analytical computation that the rational terms of all Feynman diagrams add up to zero for 
all helicity configurations. 

Our implementation of this formalism is designed for arbitrary processes with up to six 
external legs, including massive particles, and makes it possible to obtain rational terms for 
phenomenologically relevant partonic amplitudes in an automated way. Note that numer- 
ical instabilities are typically mild in terms which are free of logarithms and dilogarithms. 
The method is thus a complement to the unitarity based techniques, which lead in general 
to compact representations for coefficients of non-polynomial terms. Combining both meth- 
ods thus might be a very fruitful starting point for a highly effective method to evaluate 
complex multi-leg one-loop amplitudes. 



Acknowledgements 

First of all, we would like to thank Pierpaolo Mastrolia, because it is thanks to many 
interesting discussions with him that we focused particularly on this subject. T.B. and 
G.H. are grateful to the LAPTH Annecy for its hospitality while part of this work has 
been completed. The work of T.B. was supported by the Deutsche Forschungsgemeinschaft 
(DFG) under contract number BI 1050/1 and the Scottish Universities Physics Alliance 
(SUPA). G.H. was supported by the Swiss National Science Foundation (SNF) through 
grant no. 200021-101874. 



A. Rational parts of IR divergent integrals 

As an alternative to the approach outlined above, where IR regulated amplitudes were 
considered, one can also define form factors for the rational parts in the presence of infrared 
poles. In section 2, we first apply the operator which extracts the rational parts and then 
take the on-shell limits of the form factors F: 



lim {V + 7^) F(IR-rcgulated) 

IR L 



(A.l) 



Another possibility is to apply the operators V and TZ directly to the potentially IR- 
divergent form factors, i.e. do the operation 



{V + n) F(IR-divergent) 



(A.2) 



In an infrared finite amplitude, of course all finite terms coming from the expansion of 
e-dependent terms combined with 1/eiR poles finally have to cancel, such that the remaining 
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finite polynomial parts are identical to the ones obtained by procedure (|A.1| ), after summing 
over all contributions. In an infrared divergent amplitude, the finite remainders are related 
to the choice of the factorisation scheme. 

In order to be able to define the polynomial part of divergent amplitudes, it is necessary 
to single out the contributions which come from the expansion of the poles. The results 
of operation (A. 2) for the divergent three-point functions with Feynman parameters in the 



numerator, shown in Tables |^ to ^ below, are given in a form which allows to isolate these 
contributions immediately. 

The definition of the rational parts of expressions containing and 1/e poles is of 
course linked to the overall e-dependent factors which have been extracted. We single out 
the pole contributions in terms of ^^[13 (0, 0, X)] and U[l2]-, i.e. all double poles have been 
absorbed into the scalar three-point function with two light-like legs, depending only on 
the invariant X, and all single poles have been absorbed into I2 , where 

i2"^/2"w = ^(-xr (A.3) 



/3" = /3-(0,0,X) = (A.4) 



r(2-2e) l-2e ^ ' 

Extracting an overall factor f-p instead of rr from all integrals, we have TZ[l2] = 0, which is 
more convenient for our purposes than extracting an overall factor fr, which would imply 

nm = 2. 

For three-point functions with one non-zero invariant X, we labelled the internal prop- 
agators in such a way that 1S13 = X and S12 = ^23 = 0, where S is defined in eq. (2.4). For 



three-point functions with two non-zero invariants X and Y , we set (S23 = X and 5i3 = Y . 
Thus the integrals 13 (i, j, . . . ; 0, 0, X) are symmetric under exchange of 1 ^ 3 and the 
. . . -,0, X,Y) are symmetric under simultaneous exchange of 1 <-> 2 and X ^ Y. 
In the tables, the results for the divergent three-point functions with Feynman param- 
eters in the numerator are split up in the following way: 

mam = dsU[IS{0,0,X)]+d2U[I^]+Vi{ji}) + W{{ji}) , (A.6) 

where ds is equal to one if the integral has a double pole, and zero otherwise. The func- 
tions V{{ji}) and VF({jz}) denote the remaining finite part, where V{{ji}) is the part 
which is equal to the limit of the corresponding expression for the off-shell integral (given 
in eqs. (p. 261), (2. 27)) when one or two invariants go to zero, and W{{ji}) is the finite re- 



mainder, i.e. the difference to the on-shell limit of the corresponding off-shell integral, after 
having singled out the pole contributions. To give an example, we have 

U[I^i2, 2, 2; 0, 0, s, )] = U[I^{0, 0, si)] + ^ U[I^] + ^ 

hm U[i:^{2, 2, 2; S3, S2, si)] = - lim UiAl:^^] = - = V(2, 2, 2; si) 

19 1 
^ Wi2, 2, 2; Si) = — - V(2, 2, 2; ^i) = — . (A.7) 
9si 9si 
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The fact that W is not always zero of course does not mean that the results obtained 
by procedures (A.l) and ( [A.2| ) have to be different. The coefficients of the corresponding 
integrals as well as the number of non-zero two-point functions are different in the two 
approaches, such that after summation over all contributions making up a finite amplitude, 
the results will be the same. This has been checked explicitly by calculating the 4-photon 
and the 6-photon amplitudes in both ways. 
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V 






W[/3"] 








I^{1;0,0,X) 
IS{2;0,0,X) 
/3"(3;0,0,X) 














l3-(l,l;0,0,X) 
/3"(2,2;0,0,X) 

J3 J, U, U, v\ j 

/f(l,2;0,0,X) 
/3"(1,3;0,0,X) 
/f(2,3;0,0,X) 


1 
2X 

~^^[^2] - 2X 




1 

X 

n 
u 

1 

2X 

1 

2X 

1 
2X 





n 
u 








/|^(1,1,1;0,0,X) 
/f(2,2,2;0,0,X) 
/f(3,3,3;0,0,X) 
/f(l,l,2;0,0,X) 
7f(l,2,2;0,0,X) 
/|^(1,1,3;0,0,X) 

^3 {■^1 ^ ^ 

/|^(l,3,3;0,0,X) 
/3"(2,3,3;0,0,X) 
/3"(1,2,3;0,0,X) 


~3X^[^2] - Tlx 
W[/3"] + iiZ^[^2"] + # 
~3X^[-^2] ~ TKX 
~6X^[-^2] - 9X 

~3X^[-^2] ~ 9X 

1 

6X 

3X'^i-'2J 9X 

1 
6X 

~6X^[-^2] - 9X 

1 
6X 




2 
X 



1 

6X 

1 
2X 

1 

6X 

1 
2X 

1 

6X 

1 
6X 

1 

6X 


1 

18X 

1 

9X 
1 

18X 
1 

18X 
1 

18X 



1 

18X 


1 

18X 





Table 1: Rational and pole parts of three-point functions with two light-like legs and up to three 
Feynman parameters in the numerator. U is the operator extracting the pole and rational parts of 
the integral. The 1/e and poles have been absorbed in the terms proportional to and 
respectively. V denotes the value obtained from the corresponding expression for the off-shell 
integral, in the limit where two legs go on-shell. W is the difference h{[I] — V, where the pole terms 
have been set to zero. 
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U[I] 


V 


W 













I^{1;0,X,Y) 
IS{2;0,X,Y) 
ISiS;0,X,Y) 















/3"(1,1;0,XF) 
IS{2,2;0,X,Y) 

T^ci n X Y) 

-'3 v^) o, u, w, J J 

ISil,2;0,X,Y) 
ISil,3;0,X,Y) 
ISi2,3;0,X,Y) 


2(x-y)^ "L-'2J 1 (x-y)^ 

3y-X 7/rrnl , y 

2(x-y)2 '^H-'2J 1 (x-y)2 

n 

5^f±^ {U[I^] + 1) 
1 

^(^VL — r ; 
1 

2(y-x) 


X 

(x-y)^ 
y 

(x-y)^ 
n 

x+y 
2(x-y)^ 

1 

1 

2(y-x) 





n 






IS{1,1,1;0,X,Y) 
ISi2,2,2;0,X,Y) 
IS{S,S,S;0,X,Y) 
/f(l,l,2;0,X,F) 
/f(l,2,2;0,X,F) 
1^(1,1, 3; 0,X,F) 
/?(2,2,3;0,X,y) 
73"(l,3,3;0,X,y) 
73"(2,3,3;0,X,y) 
/3"(1,2,3;0,X,F) 


iix2-7xy+2y2 7yrrni i 37x2-8xy+y2 


x(6X-y) 


1 


6{x-y)i* '^'L-'2J"r i8(x-y)^ 
iiy2-7xy+2x2 7/rrni i 37y2-8xy+x2 


3(X-y)a 

y(6y-x) 


i8(x-y) 
1 


6(y-x)3 ^-'1^2 J 1 i8(y-x)a 


-2x2-5xy+y2 7/[rni 1 -5x2-iixy+y2 


3(y-X)a 



-3x2-8xy+y2 


i8(y-x) 


1 


6(x-y)^ ^n^2\ 1 9(x^y):* 
-2y2-5xy+x2 7/frni 1 -5y2-iixy+x2 


6(x-y)i' 
-3y2-8xy+x2 


i8(x-y) 
1 


6(y-x)^ ^n^2\ 1 9(y-x)a 

-3X+y 
6(x-y)2 

-3y+x 
6(x-y)^ 

1 

6(x-y) 
1 

6(y-x) 

(x+y) 
6(x-y)^ 


6{y-x)i' 

-3X+y 
6(x-y)2 

-3y+x 
6(x-y)^ 

1 

6(x-y) 

1 

6(y-x) 

(X+y) 
6(x-y)^ 


i8(y-x) 












Table 2: Rational and pole parts of tlircc-point functions with one light- like log and up to three 
Feynman parameters in the numerator. U is the operator extracting the pole and rational parts 
of the integral. V denotes the value obtained from the corresponding expression for the off-shell 
integral, in the limit where one leg goes on-shell. W is the difference U[I]—V where the pole terms, 
absorbed into have been set to zero. 
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^J\T^ 


V 


W 


I^+\0,0,X) 




-Mm - \ 





/3"+2(l;0,0,X) 

/3"+2(3;0,0,X) 


-\m] - \ 

-Wim - 1 


-^\m - IE 
-M^2\ - \ 







I^+\0,X,Y) 


-Mm - h 


-Mm-h 





ltHl-0,X,Y) 
I^+\2;0,X,Y) 


qI^[12\^ 18(X-y) 
gM[i2j -1- I8(y-X) 


6^[^2\^ 18{X-Y) 
gM[i2j -1- i8(X-y) 

-WE] -I 








Table 3: Rational and pole parts of (n + 2)-dimensional three-point functions with one or two 

light-like legs and up to one Feynman parameter in the numerator. U is the operator extracting 
the pole and rational parts of the integral. Note that the poles in T^^^ are of ultraviolet nature, 
but we do not distinguish the nature of the poles denoted by 'P[l2]- V denotes the value obtained 
from the corresponding expression for the ofF-shcU integral, in the limit where one or two legs go 
on-shell. W is the difference h{[I] — V, where the pole terms, absorbed into W[/2 ], are set to zero. 
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